A Compound Josephson Junction Coupler for Flux Qubits With Minimal Crosstalk 



0^ 
O 
O 



R. Harris/' * T. Lanting/ A.J. Berkley,^ J. Johansson,^ M.W. Johnson/ 

P. Bunyk/ E. Ladizinsky/ N. Ladizinsky/ T. Oh/ and S. Han^ 

^D-Wave Systems Inc., 100-4401 Still Creek Dr., Burnaby, BC V5C 6G9, Canada) 
^Department of Physics and Astronomy, University of Kansas, Lawrence KS, USA 

(Dated: July 17, 2009) 

An improved tunable coupling element for building networks of coupled rf-SQUID flux qubits has 
been experimentally demonstrated. This new form of coupler, based upon the compound Joseph- 
son junction rf-SQUID, provides a sign and magnitude tunable mutual inductance between qubits 
with minimal nonlinear crosstalk from the coupler tuning parameter into the qubits. Quantita- 
tive agreement is shown between an effective one-dimensional model of the coupler's potential and 
measurements of the coupler persistent current and susceptibility. 



o 
o 

Oh: 



o 



> 
oo 
m 

o 

On 
O 



X 



PACS numbers: 85.25.Dq, 03.67.Lx 

The choice of architecture of a prototype solid state 
quantum information processor is primarily driven by 
the algorithm that the designer wishes to implement. 
Within the field of superconducting quantum devices, at 
least two distinct architectures have arisen. Gate model 
algorithms require qubits with long-lived excited states 
and dynamic couplings. Recent efforts have focused upon 
chargc-likc [1] and phase [2] qubits coupled to microwave 
resonators. Adiabatic quantum algorithms [3, 4] require 
qubits whose groundstate encode binary variables and 
static couplings. One implementation involves a network 
of inductively coupled flux qubits [5] . The Hamiltonian 
for this architecture is that of a quantum Ising spin glass, 
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where = 2\If\ $f and arc the bias and tunneling 
energy of qubit i, respectively, and = My is 
the coupling energy between qubits i and j. Here, \If\ 
represents the magnitude of the qubit persistent current, 
$^ is an external flux bias and Mij is a mutual induc- 
tance. A programmable processor would require in-situ 
tunable <i>f and Mij. Inductive coupling could also be 
useful in other quantum computation schemes in which 
the flux qubit 's persistent current basis is nearly con- 
current with the computation basis [6]. More involved 
parametric coupling schemes are needed if the flux qubits 
are biased to their optimal points where the energy and 
persistent current bases arc orthogonal [7]. 

The authors of Ref. [8] proposed the use of an rf- 
SQUID to implement tunable Mij. Experiments on sys- 
tems of coupled flux qubits verified that such couplers 
did perform as anticipated [9]. However, additional work 
not reported in the literature revealed two serious defi- 
ciencies: First, the tuning mechanism involves threading 
flux through the rf-SQUID loop, thus inducing a large 
persistent current Ip that, in turn, biases the qubits. 
This is a significant problem if the qubit biases need to 
be controlled to high precision atop what can be a very 




FIG. 1: (Color online) Schematic of a CJJ rf-SQUID. 



large nonlinear crosstalk imparted by the coupler. Sec- 
ond, Alij = can only be achieved if /3 = 2ttLIc/^o < 1, 
where L and Ic are the rf-SQUID inductance and criti- 
cal current, respectively. On the other hand, in order to 
achieve appreciable non-zero coupling it proved necessary 
to design devices with (3 > 0.9. Such devices were acutely 
sensitive to fabrication variations, where higher than ex- 
pected Ic could make My = unattainable. Thus the 
rf-SQUID coupler proved troublesome in practice. Note 
that the dc-SQUID [8, 10] and the unipolar [11] couplers 
suffer from similar deficiencies. The challenge was then 
to design a tunable My that invokes minimal Ip and is 
robust against fabrication variations. 

Our efforts to identify a satisfactory coupler design led 
us to consider the compound Josephson junction (CJJ) 
rf-SQUID, as depicted in Fig. 1. The CJJ rf-SQUID com- 
prises a superconducting loop of inductance L that is in- 
terrupted by a smaller loop of inductance Leo containing 
two Josephson junctions with critical currents Id and Ic2- 
Devices with this general topology can take on a num- 
ber of guises: The CJJ rf-SQUID can be used as a qubit 
when designed with relatively low Ic = Id + Ic2 ^ 3 fiA, 
L chosen such that /3 < 2, low net capacitance across 
the junctions C and biased such that its potential en- 
ergy is bistable [12, 13]. The CJJ rf-SQUID can also be 
used as a latching readout when designed with a sub- 
stantial Lc > 10 ^A, L chosen such that /3 > 1.2, large C 
and having its potential swept from being monostable to 
bistable [14]. We focus herein upon a CJJ rf-SQUID de- 
signed with a modest Lc, L chosen such that (3 < 1.2 
and operated with a monostable potential. In Fig. 1 
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the closed loops are subjected to external flux biases 
*act = *o'Pact/27r and 3>^„ = <S>oVco/'^^- The symbols 
have been chosen to indicate that ^^^^ is an actuator for 
perturbing the device and represents the control sig- 
nal. Let the phase across the junctions be ipi and ip2- 
The Hamiltonian for this device can be written as 
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FIG. 2: (Color online) Schematic of a CJJ rf-SQUID coupler 
interacting with two CJJ rf-SQUID qubits. 



where Ci and Eji = Id^o/^n represent the capacitance 
and Josephson energy of junction i, respectively, and 
[$o<y5j/27r, Qj] = ihSij. The inductive terms originate 

from the two closed loops with n £ {co, act}, L^ct = which can be solved for (pact given arbitrary i^^cti thus 
L+Lco/4: and J7„ = ($o/27r)^/i„. The actuator and con- 
trol loop phases are defined as (^act = {^i + ^2) /2 and 
ifco = ^1 — ^2, respectively. Hamiltonian (2) can be re- 
duced to an effective 1-dimensional system if L^ct ^ Leo 
because the plasma energy of the control loop will then 
be much higher than that of the actuator loop. Setting 
fco = fco ^^'^ combining the Josephson terms, 



yielding /fct(*^act' *^'co)- Differentiating Eqs. 5 and 6 with 
respect to ^^^.^ then yields x^^^: 
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where Ic± = Id ± Ic2 and Cp = Ci -|- C2 . Hamiltonian 
(3) is homologous to that of an rf-SQUID whose single 
junction possesses a critical current that is a function of 
(fi^^ and whose phase has been shifted by if'^^^- 

Let the device described by Eq. (3) be connected to 
two qubits via mutual inductances Mco,i and A/co,2- The 
mutual inductance between the qubits will be 
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where x^^-* = dllctl '^'^tct represents the first order (lin- 
ear) susceptibility of the coupler [8] and the persistent 
current fiowing about the coupler actuator loop is 
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If V{ipact) is monostable and the first excited state can 
be neglected, then one can replace the operator (pact by 
the value for which y is a minimum [dV/dLp^ct = 0): 



Vs.ct - vlct + /3cff sin ((^a 
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Equation (7) is similar to Eq. (10) of Rcf. [8], albeit PcS 
is a function of ip^^ and junction asymmetry results in a 
<p^Q-dependent phase shift in the cosine terms. 

While rf-SQUID and CJJ rf-SQUID couplers possess 
similar expressions for x'-^^ the latter holds two advan- 
tages: First, the CJJ coupler can be operated with <I>act = 
and tuned via $^0- If Ic-llc+ < 1, then (^^^.j < 1 and 
Eq. (6) yields (ySact ~ 0. Equation (5) then predicts that 
/act ~ 0- Thus the CJJ coupler need not invoke large per- 
sistent currents (on the order of Ic+) when being tuned. 
Second, the CJJ coupler is usable over the range of ^'^^ 
for which -min [1, /3eff(0)] < (iesivlo) < /5eff(0) when 
'^Ict ^ 0' where the lower bound has been imposed by the 
condition that I^(v?act) be monostable. Thus the utility 
of the CJJ coupler is not compromised if /3eff(0) > 1. As 
such, this device is robust against fabrication variations. 

To test the CJJ rf-SQUID coupler, we fabricated a 
circuit containing 8 CJJ rf-SQUID flux qubits [12, 13], 
each inductively coupled to its own hysteric dc-SQUID 
readout [15], and connected by a network of 16 CJJ rf- 
SQUID couplers. The chip was fabricated from an oxi- 
dized Si wafer with Nb/Al/Al203/Nb trilayer junctions, 
four Nb wiring layers capped with SiN and separated by 
planarized PECVD Si02. The chip was mounted to the 
mixing chamber of a dilution refrigerator regulated at 
T = 40 mK inside a Sn superconducting magnetic shield 
with a residual field in the vicinity of the chip < InT. 
External current biases were provided by room tempera- 
ture current DACs whose outputs were low pass filtered 
with /c ~ 5 MHz using a combination of lumped element 
and copper powder filters secured to the mixing chamber. 

A schematic of a coupler and a pair of qubits is de- 
picted in Fig. 2. The coupler is controlled via bias cur- 
rents that are coupled to the device through mutual in- 
ductances Mco and Mact , respectively. These give rise to 
the fluxes and 'I'acf The qubits are controlled via 
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FIG. 3; (Color online) a) Single qubit measurement wave- 
forms, b) Sequential annealing waveforms. Source (detector) 
qubit waveforms denoted as dashed (solid) lines. 



fluxes ?>^jj„ and 



(a = 1,2) as described in Ref. [13]. 
The qubits interact with the coupler via mutual induc- 
tances Mco,a- For brevity, we present results from a sin- 
gle coupler in this paper and note that Afoff(*&co) 
identical to < 5% for all 16 couplers on this chip. For the 
particular coupler described herein, the relevant qubit 
critical currents were /^^ = 3.25 ± 0.01 and qubit in- 



ductances were L, 



,1(2) = 290(308) ± 5 pH when = 0. 
The flux waveforms used to obtain Mact are depicted 
in Fig. 3a. In this case, was held constant while the 
detector qubit (a = d) was annealed in the presence of a 
pulse on <i>act (*) of amplitude $act and a pulse on <^^^ {t) 
of amplitude The sequence involves initializing the 
qubit in a monostable potential with no net flux biases 
(i), setting ^^^.^ and (^^^ (ii), raising the detector qubit 's 
tunnel barrier to maximum height (iii) and then return- 
ing and <I>^jj to zero prior to reading the state of the 
detector qubit (iv). The result of this process is that the 
state of the detector qubit becomes trapped in one of its 
two counter-circulating persistent current states, denoted 
as |0) and |1). Repeating this sequence to gather statis- 
tics then yielded the probability of finding the detector 
in |0), Pq- Using software feedback, we adjusted <i>^ to 
track the displacement of the detector's degeneracy point 
^>°Q,, defined as the bias for which Pq = 1/2, to within 
±0.02 m$o- We have defined <i>°^ = with respect to the 
degeneracy point obtained with — and <I>act — 0- 
Mapping <I>°^ versus the current bias driving 4>^^^ yielded 
a modulation with period A/act, from which we obtained 

Mact = $o/A/act = 1-77 ± 0.01 pH. 

To obtain Mco we again used the flux waveform pat- 
tern depicted in Fig. 3a but with $act toggled between 
±5m$o- Taking the difference in (f>° ^ between the two 
polarizations, we tracked the amount of coupled flux 
= 2Mco,ax'"^\^to)'^\ct versus the bias driving <^l^. 
The results yielded a period A/co, from which we ob- 
tained Afco = *o/A/co = 3.43 ± 0.03 pH. 

With the coupler biases calibrated, we proceeded with 
measuring Adcsi^^o)- To do so, we used the 2-qubit flux 



bias sequence depicted in Fig. 3b in which one qubit 
served as a flux source {a = s) and the other acted as 
a flux detector [a = d). This process, referred to as 
sequential annealing, involved initializing both qubits in 
monostable potentials with no net flux biases (i), set- 
ting = = ±5m$o and raising ^cjj.s t''^? 
source qubit in either |0) or |1) (ii), using software feed- 
back to adjust <i>^ (iii) and raising the detector qubit's 
tunnel barrier to trap its state (iv). Finally, both 
were returned to zero prior to reading the state of the 
detector qubit (v). The relative change in <I>°^ between 
the two polarizations of the source qubit then yielded the 
flux X^"^ = 2Moff($^o) For each qubit, could be 
directly inferred from measurements obtained with its 
readout. Measurements of Mcs{^co) ^^^^ shown in Fig. 4. 
Using any three of {x'"^\ x''^\ x'^^ ) , one could also 
solve for Mco,i = Mco,2 = 17.5 ± 0.2 pH and x^^K^co)- 

A key motivation for developing the CJ J coupler was to 
minimize the impact of coupler settings upon qubit prop- 
erties. Plots of the relative change in qubit degeneracy 
point (f>°|_, = Mcca/act vcrsus are shown in Fig. 5a. 
The qubit inductance Lqa will also be a function of ^^o- 
Let the change in inductance be defined as 
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We have quantified this effect by measuring | | versus 



<f>Jo. For ^cii.Q ^ and (f>^^ ~ one can use an 



approximation, similar to that used to obtain Eq. (5), to 
write a pair of expressions for 
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FIG. 4: (Color online) CJJ rf-SQUID coupler effective mutual 
inductance versus control flux. Solid (hollow) points corre- 
spond to Of = 1(2) acting as flux source and a = 2(1) as flux 
detector. Solid curve is from a simultaneous best fit of these 
data to Eqs. (4) and (7) and those in Fig. 5a to Eq. (5). 
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FIG. 5: (Color online) Deviation of qubit parameters as a 
function of coupler control flux: a) Qubit degeneracy point. 
Solid curve is from a simultaneous best fit of these data to 
Eq. 5 and those in Fig. 4 to Eq. (7). b) Qubit inductance. 
Solid curve is prediction using the best fit parameters. 



Given the calibrated values of for each qubit, it was 
then possible to extract Lgai^^a) from |/^(<I>co)l- The 
resultant 6Lqa = Lqa{^^^) ~ Lqa{0) are shown in Fig. 5b. 

The data shown in Figs. 4 and 5a have been simul- 
taneously fit to Mco,iMco,2X^^^ [Eqs. (4) and (7)] and 
Mco,aIact [Eq- (5)], respectively. In order to yield a high 
quality fit. we allowed for small flux offsets in both cou- 
pler loops and a small linear crosstalk from the con- 
trol bias to the actuator loop: — and 



act 



7$^Q — $acf The best fit was obtained with 



Ic+ = 8.4±0.3M, Ic-/Ic+ = (6.7 ± 0.9) x 10"^ 
iact = 88±2pH, = 16±lm$o, *act = 3±lm$o and 
7 = (6 ± 1) X lO"''. Note that direct coupling between 
qubits would add a positive constant to A/cff [Eq. (4)]. 
A key consequence would be that the difference between 
the two values of where M^s = would be greater 
than <I>o, as predicted by Eq. (7). The best fit curve in 
Fig. 4 suggests that A^^^Im^^^^q w 0.998 $o, therefore 
the direct coupling appears to be negligible. 

The solid curve in Fig. 5b represents the predicted 
dLqa using the best fit parameters. Given the agree- 



ment between theory and experiment, one can conclude 
that Eqs. (3)^(7) correctly model the physics of this 
device. Note that over the bias range shown that the 
qubit degeneracy points shift by ~ 2m$o *C Mco,aIc+ ~ 
70 m$o- Consequently, one can conclude that the nonlin- 
ear crosstalk from coupler to qubit is substantially less 
than that encountered while tuning a comparable sin- 
gle junction rf-SQUID coupler [8]. According to Eq. 5, 
this undesirable effect could be reduced to negligible lev- 
els by improvements in fabrication uniformity (smaller 
Ic-/Ic+)- Achieving lower 7 in future designs will also 
help realize further reductions in nonlinear flux offsets. 

Changes in Lqa are of consequence if the properties of 
multiple qubits need to be synchronized to high preci- 
sion [13]. Custom tuned qubit CJJ flux offsets provide 
one means of mitigating this undesirable effect [13]. Al- 
ternate qubit designs which contain an in-situ tunable 
inductance for ballast constitute a second solution. 

Conclusions: A compound Joscphson junction rf- 
SQUID coupler suitable for building networks of coupled 
flux qubits has been demonstrated. This coupler provides 
both sign and magnitude tunable mutual inductance in 
a manner that invokes minimal nonlinear crosstalk from 
the coupler tuning parameter to the qubits. Furthermore, 
this crosstalk can be reduced to negligible levels with im- 
proved fabrication uniformity and subtle improvements 
in device layout. Modulation of the qubit inductance 
via changes in the coupler settings has been character- 
ized and shown to be predictable using an effective one- 
dimensional model of the coupler potential. 
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